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Motivated by the recent experimental evidence of commensurate surface charge-density-waves 
(CDW) in Pb/Ge(lll) and Sn/Ge(lll) y^-adlayer structures, as well as by the insulating states 
found on K/Si(lll):B and SiC(OOOl), we have investigated the role of electron-electron interactions, 
and also of electron-phonon coupling, on the narrow surface state band originating from the outer 
dangling bond orbitals of the surface. We model the -^3 dangling bond lattice by an extended 
two-dimensional Hubbard model at half-filling on a triangular lattice. The hopping integrals are 
calculated by fitting first-principle results for the surface band. We include an on-site Hubbard 
repulsion U and a nearest-neighbor Coulomb interaction V, plus a long-ranged Coulomb tail. The 
electron-phonon interaction is treated in the deformation potential approximation. We have explored 
the phase diagram of this model including the possibility of commensurate 3x3 phases, using 
mainly the Hartree-Fock approximation. For U larger than the bandwidth we find a non-collinear 
antiferromagnetic SDW insulator, possibly corresponding to the situation on the SiC and K/Si 
surfaces. For U comparable or smaller, a rich phase diagram arises, with several phases involving 
combinations of charge and spin-density- waves (SDW), with or without a net magnetization. We 
find that insulating, or partly metallic 3x3 CDW phases can be stabilized by two different physical 
mechanisms. One is the inter-site repulsion V, that together with electron-phonon coupling can 
lower the energy of a charge modulation. The other is a novel magnetically-induced Fermi surface 
nesting, stabilizing a net cell magnetization of 1/3, plus a collinear SDW, plus an associated weak 
CDW. Comparison with available experimental evidence, and also with first-principle calculations 
is made. 

I. INTRODUCTION 

Pb and Sn \/3 x \/3 adlayer structures on the (111) surface of Ge have recently revealed a reversible charge density 
wave (CDW) transition to a low temperature reconstructed 3x3 phase.ErQ A half-filled surface state band makes 
the high temperature phase metallic. The low temperature phase is either metallic - as seems to be the case for 
Sn/Ge(lll) - or weakly gapped, or pseudo-gapped, as suggested for Pb/Ge(lll). 

Related isoelectronic systems, like the v^-adlayer of Si on the (0001) surface of SiCo and on K/Si(lll):B,El show 
a clear insulating behavior, with a large gap, no structural anomalies, no CDWs, and no transitions, at least to our 
present knowledge. 

These adsorbate surfaces are altogether mysterious. The very existence of a \/?> adsorbate phase, with coverage 
1/3, is puzzling. For isoelectronic Si on Si(lll), or Ge on Ge(lll), for instance, there exists no such phase. The stable 
low-coverage phase are 7x7 and c(2 x 8) respectively, whose coverage is instead close to 1/4. They are made up of 
2x2 basic building blocks, each with one adatom saturating three out of four first-layer atoms, and one unsaturated 
first-layer atom, the "restatom" . In this adatom-restatom block, the nominally unpaired electron of the adatom and 
that of the restatom pair off together, giving rise to a stable, fully saturated, insulating surface. By contrast, the 
•\/3 x -\/3 phases with 1/3 coverage are very common for trivalent adsorbates, such as Ga and In, and for pentavalent 
ones like As, on the same (111) surfaces. These adatoms lack the unpaired electron, and can therefore lead to a fully 
saturated insulating surface without the need for any restatoms. 

A \/3 x \/3 adsorbate phase of tetravalent adatoms is bound by construction to possess one unpaired electron per 
adatom, giving rise to a very destabilizing half-filled metallic surface state band. Seen in this crude light, it is a puzzle 
why this kind of coverage should constitute even only a locally stable state of the surface. 

Lookingjuacfi closely, we may speculate that SiC(0001)El and K/Si(lll):B,Q most likely Mott-Hubbard 
insulatorsJiy~t3ia are perhaps "stabilized" by Coulomb repulsions, so large to make it anyway difficult for electrons to 
move. For the more innocent-looking, less correlated, Pb/Ge(lll) and Sn/Ge(lll), this argument is less obvious, and 
the puzzle remains. The function of the 3x3 CDW state - whatever its real nature - most likely serves the function 
of stabilizing these otherwise unstable surfaces at low temperatures. Nonetheless, the CDW periodicity chosen by 
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the surface CDW -3x3, meaning a super-cell of adatoms - is not at all evident. In fact, it replaces a 

supposedly unstable state characterized by an odd number of electrons/cell (three), with another where the electron 
number (nine) is, alas, odd again. 

Be all that as it may, there is little doubt that the main factor driving the phenomena on all these surfaces, appears 
to be precisely the half-filled- and extremely narrow - surface state band. We thus begin with a discussion that in 
principle encompasses all the y/3 x \/3 tetravalent adsorbed surfaces. 

We believe the following points to be of general validity: 

i) Poor nesting. . 

Two-dimensional Fermi surface (FS) nesting in the half-filled sou-face statestLS has been repeatedly irnpked as the 
driving mechanism for the CDW instability in^ha-pase of Pb/GeJj'El but excluded for the case of Sj^/Xic.ErEj However, 
by inspecting either photoemission k(E) datajjuu and existing first-principle (LDA) calculationsEniStil of the surface 
half- filled band (the "adatom dangling bond band" ) , we fail to detect a particularly good nesting of the two-dimensional 
FS at the surface Brillouin zone (BZ) corner K = (47r/3a,0). The wavevector-dependent susceptibility generated by 
the calculated band structure, in particular, has no especially large value at this k-point, and rather peaks elsewhere 
(see inset in Fig. 0). To be sure, there is nothing preventing in general a good nesting at K = (47r/3a,0), or any 
other k-point. However, insofar as the surface state band is really lying in a bulk gap at each single k-point, it should 
be with good accuracy - by simple state counting and charge neutrality - precisely half filled. This implies that the 
filled and empty state areas should be equal. Hypothetical Fermi surfaces with this kind of shape and good nesting 
at K = (47r/3a,0) do not appear to be compatible with an integer electron number. We thus believe lack of perfect 
nestingto be the case for both Pb/Ge as for Sn/Ge. ._. 

Fig. hi showing a tight binding fit to the LDA surface band dispersion for the test-case of Si(lll)/Siji3 as well as 
the corresponding FS and Lindhard density response function Xo(q), 

, v f dk n k - n k+q 

JBZ ( 27T ) £ k+q - £k 

rik and ek being the occupation number and energy of an electron with Bloch momentum k, provides a concrete 
illustration of these statements. We note, in passing, that a strong nesting at K is, on the contrary, automatically 
guaranteed if the surface band acquires a tuaiform magnetization in such a way that the densities of up and down 
electrons become, respectively, 2/3 and 1/3.E3 The majority spins would then fill the region external to the reduced 
BZ in Fig. |], and their FS would be strongly nested. This suggestion, which turns out to be correct at the mean-field 
level, points into the direction of a possible role played by magnetism in these systems. 

ii) Importance of electron-electron interactions. 

The width W of the surface band is relatively small: W « 0.5 eV for Pb and Sn/Ge(lll), W » 0.3 eV for SiC(OOOl). 
Moreover, this band is half-filled. These facts call for a careful consideration of electron-electron interactions, as well 
as of electron-phonon (e-ph), as possible sources of instability. The importance of electron-electron interaction is un- 
derlined by the different phenomenology of SiC(OOOl) and K/Si(lll):B with respect to Pb-Sn/Ge(lll). The stronger 
insulating character of the former surfaces parallels closely their stronger electron-electron repulsions, connected both 
with more localized surface Wannier functions (see later on), and with reduced screening, due to larger bulk semicon- 
ducting gaps. 

Hi) Weakness of LDA calculations for ground state prediction. 

LDA electronic structure calculations - an extremely well tested tool in many areas- are certainly suitable for a 
weakly interacting system, such as the bulk semiconductor, or a passivated semiconductor surface. They are less 
reliable, especially when they do not include spin, in predicting the stable state and the instabilities of a narrow band 
system. For instance, the phenomenology of SiC(OOOl) - suggesting a Mott-Hubbard insulator.^, is unreproducible 
by LDA. The onset of a CDW on Sn/Ge(lll)is also not predicted by recent LDA calculationsBEJ While there is no 
reason to doubt the basic credibility of the one-electron band energies obtained from these Kohn-Sham equations, the 
mean-field treatment of interactions, the screened local exchange, and especially the neglect of magnetic correlations 
are the standard source of problems with LDA. As a consequence, it will be necessary to worry more substantially 
about interactions, and to use methods which, even if mean-field, permit the inclusion of strong correlations, including 
magnetic effects. 

iv) Interaction-driven mechanisms for 3x3 CDW instabilities. 

There are several different couplings which the surface electrons, as they hop weakly between a surface adatom site 
and another, experience and can influence the formation of the CDW, or of an insulating ground state: a) on-site, 
and nearest-neighbor (n.n.) inter-site electron-electron repulsion; b) on-site effective attraction (negative Hubbard-C/ 
term) of electron-phonon origin. 

Because of poor nesting, electron-phonon alone is unlikely to drive the 3x3 CDW. At weak coupling, the suscepti- 
bility peak in Fig. [I] would rather drive an incommensurate periodicity. At strong coupling, the frustration associated 
to the triangular lattice, will favor, in general, a superconducting ground state over a CDW phase (see Appendix) £3 
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On the other hand, the electron-electron interaction, both on-site and, independently, nearest neighbor, naturally 
suggests, as we shall see later, the 3x3 surface periodicity, which is found experimentally. 

The approach we will take is based on an extended Hubbard-Holstein model. It is by necessity a "non-first- 
principle" approach and, as such, has no strong predictive power. However, it is made more realistic by using 
parameters extracted from first-principle calculations, and we find it very helpful in clarifying the possible scenarios 
as a function of the strength of electron-electron interactions. Because of this rather qualitative use, we will make no 
attempt to push the accuracy of treatment of this model to a very high level of sophistication. The basic tool will be 
the unrestricted Hartree-Fock approximation. Although mean field, it allows magnetic solutions, favored by exchange 
which is unscreened. 

II. MODEL 

Each tetravalent adatom on a (111) semiconductor surface carries a dangling bond - an unpaired electron in an 
unsaturated orbital. In the V3 x a/3 structure, the dangling bonds of the adatoms give rise to a band of surface states 
which lies in the bulk semiconductor gap.EHLj By electron counting, such a band is half-filled. Our basic starting 
is the quantitatively accurate surface state band dispersion et which one calculates in gradient-corrected LDA 
It is shown in Fig. [j] for the case of Si/Si(lll). The solid and dashed lines in Fig. | are tight-binding fits to the 
LDA results obtained by including, respectively, up to the 6 th and up to the 2 nd shell of neighbors. The fit with 
hopping integrals t\,t2, ■ ■ ■ ,tg is quite good. Less good, but qualitatively acceptable, is the fit obtained using only 
nearest neighbor (n.n.) and next-nearest neighbor (n.n.n.) hopping integrals t\ and £2- The Fermi surface (FS) for 
the half-filled surface band is shown in the upper inset of Fig. [y. It is important to stress that the FS does not show 
good nesting properties at the wavevector q = K (the BZ corner). This feature is shared by all LDA calculations on 
similar systems BEiffi Albeit small, the bandwidth W of the surface band is much greater than one would predict by 
a direct overlap of adatom dangling bonds, as the adatoms are very widely apart, for instance about 1A on Ge(lll). 
Hopping is indirect, and takes place from the adatom to the first-layer atoms underneath, from that to a second-layer 
atom, then again to a first-layer atom underneath the other adatom, and from there finally to other adatom dangling 
bond. Thus, when expressed in terms of elementary hopping processes between hybrid orbitals, electron hopping 
between two neighboring adatom dangling bonds is fifth order. As a result, the final dispersion of the surface state 
band strongly parallels that of the closest bulk band, the valence band. Correspondingly, hybridization effects of the 
dangling bond orbitals with first, second, and even third, bulk layer orbitals are strong, as shown by the extension 
into the bulk of the Wannier orbital associated to the LDA surface band (Fig. |^) . 

In spite of this, we can still associate to every adatom a Wannier orbital and write the effective Hamiltonian for 
the surface band as follows: 

BZ 

H = J2J2 ^claC Ka + H ph + ife-ph + Hint , (1) 
k a- 

where is the Fourier transform of the Wannier orbital, namely the surface state in a Bloch picture. The sum 
over the wavevectors runs over the surface BZ. H mi includes correlation effects which are not correctly accounted for 
within LDA, which we parametrize as follows: 

Hi n t = ^y^n r ,i-n r ,4 + i ^2 Vr-r'(n T - l)(»v - 1) ■ (2) 

r r^r' 

Here U is an effective repulsion (Hubbard-[/) for two electrons on the same adatom Wannier orbital, and V T - r > is the 
direct Coulomb interaction between different sites r and r'.EZl Let V be the n.n. value of 14- r', which is, clearly, the 
largest term. We have considered two models for V r - r r. a model (A) in which we truncate V r - r i to n.n., and a model 
(B) in which V r -r' has a long range Coulombic tail of the form 

V , 



where a is the n.n. distance. The results for model B are qualitatively similar to those of A, and will be only briefly 
discussed later on. In other words, even if most of the detailed results in this paper will be base on the n.n. V r _ r ', 
their validity is more general. 

LDA estimates of the bare coulomb repulsion U a and V a between two electrons respectively on the same and 
on neighboring Wannier orbitals are - for our test case of Si(lll) /Si - of about 3.6 eV and 1.8 eV respectively.!^ 
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Screening effects by the the underlying bulk are expected to reduce very substantially these repulsive energies. An 
order of magnitude estimate for U and V is obtained by dividing their bare values by the image-charge screening 
factor, (e + l)/2 « 6, yielding, for Si, U — 0.6 eV (10ti), and V-.— 0.3 eV (5ii). Corresponding values would be 
somewhat smaller for Ge(lll), in view of a very similar dispersiorO and of a ratio of about 4/3 between the dielectric 
constants of Ge and Si. SiC(OOOl), the opposite is true. The surface state band is extremely narrow, of order 0.3 
eVt2, while the bulk dielectric constant is only about 6.5. 

As for the e-ph interaction, in principle both the on-site Wannier state energy and the hopping matrix elements 
between neighbors depend on the positions of the adatoms. Within the deformation potential approximation, we 
consider only a linear dependence of the on-site energy from a single ionic coordinate (for instance, the height z r of 
the adatom measured from the equilibrium position), and take 

#c- P h = -922 Zr ( Ur ~ ^ ' ( 3 ) 

r 

with g of the order of w 1 eV/A The free-phonon term will have the usual form 

^=2^(^ + 5) > ( 4 ) 

where is the phonon annihilation operator, and Hlo^ a typical phonon frequency of the system, which we take to 
be about 30 meV, independent of k. 



III. PHASE DIAGRAM: SOME LIMITING CASES 



Preliminary to the full treatment of Sect. IV, we consider first the purely electronic problem in the absence of 
e-ph interaction. We start the discussion from particular limiting cases for which well-controlled statements, or at 
least intuitively clear ones, can be made, without the need of any new specific calculations. In the Appendix we will 
also consider, because it is useful in connection with the electron-phonon case, the unphysical limit of strong on-site 
attraction (large and negative U). 



A. Large positive U : the Mott insulator. 

For U ^> V, W, the system is deep inside the Mott insulating regimeE! The charge degrees of freedom are frozen, 
with a gap of order U. The only dynamics is in the spin degrees of freedom. Within the large manifold of spin 
degenerate states with exactly one electron per site, the kinetic energy generates, in second order perturbation theory, 
a Heisenberg spin-1/2 antiferromagnetic effective Hamiltonian governing the spin degrees of freedom, 

Heft — ^2 ' ' 

(ij) 

with J i:j = 4|^| 2 /C/.@ 

For our test case of Si(lll) /Si, the values of the hoppings are such that J\ as 20 meV, J2/J1 ~ 0.12 while 
the remaining couplings J3, • • • are very small. Antiferromagnetism is frustrated on the triangular lattice. Zero 
temperature long range order (LRO) - if present - should be of the three-sublattice 120°-Neel type, which can be also 
seen as a commensurate spiral spin density wave (s-SDW). 

Because it does not imbalance charge, this state is not further affected by electron-phonon coupling. 

In summary, we expect for large values of U a wide-gap Mott insulator with a s-SDW (spins lying in a plane, 
forming 120° angles), a 3 x 3 magnetic unit cell, but uniform charge (no-CTJW). This is, most likely, the state to be 
found on the Si-terminated and C-terminated SiC(0001) surface at T=0EIEi 



B. Strong inter-site repulsion: an asymmetric CDW with three inequivalent sites. 

The e-ph coupling can effectively reduce U, but not V. Therefore, it is of interest to consider the hypothetical 
regime W < U <C V. When the first-neighbor electron-electron repulsion V is large the system, in order to minimize 
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the interaction energy, will prefer a 3 x 3 CDW-like ground state, with two electrons on one sublattice (A), a single 
electron on another sublattice (B), and zero electrons on the third sublattice (C) (see Fig. ||). These states are still 
highly degenerate (in the absence of hopping) due to spin degeneracy for the single unpaired electron on sublattice 
B. A gap U separates these states from the lowest-energy excited configurations (see Fig. |J). The spin degeneracy 
can be removed in second-order perturbation theory, owing to t% , which leads to an effective spin- 1/2 Heisenberg 
Hamiltonian within sublattice B, 

sublattice B 

HeS = J S '< ' S r 3 > ( 6 ) 

m 

with a weak antiferromagnetic exchange constant J — 4^/^-0 Summarizing, we expect in this regime a strong 3x3 
asymmetric CDW (a-CDW) with three inequivalent sites {<p p « 7r/6, see below), and a spiral 3\/3 x 3\/3 SDW, 
governing the unpaired electron spins, superimposed on it. Notice that, while the charge periodicity is 3 x 3, the 
actual unit cell is larger, i.e., 3-\/3 x 3^/3. Despite having the correct charge periodicity, namely 3x3, this a-CDW 
is not compatible with the experimental findings on Pb-Sn/Ge, which is a symmetric CDW. We conclude that the 
low-temperature CDW state of these systems is not completely dominated by V. 



IV. MEAN-FIELD THEORY. 



In order to get a more complete picture of additional phases for smaller U, and of the possible phase diagram 
of the model we now turn to a quantitative mean field theory analysis. The first issue is to include the possibility 
of magnetic correlations. For small values of the interactions U and V, the Stoner criterion can be used to study 
the possible magnetic instabilities of the paramagnetic metal obtamcd from LDA calculations. The charge and spin 
susceptibilities are given, within the random phase approximation,!^] by 

2Xo(q) 



xc(q) = 



l + ([/ + 2V q )Xo(q) 



where Xo is the non-interacting susceptibility per spin projection, and both factors of 2 account for spin degeneracy. 
The divergence of xs is governed, in this approximation, by U only. Since Xo(q) is finite everywhere, a finite U is needed 
in order to destabilize the paramagnetic metal. The wavevector q* at which xs first diverges, by increasing U, is in 
generalincommensurate with the underlying unit cell. The instability is towards an incommensurate, metallic, spiral 
SDW.El Fig. [I] shows that, in our case, q* = (1.32AT,0) (with K = 4n/3a, the BZ corner). We get Uf F jt x rs 3.7. 
(The other maxinmm of Xo at q = (0.525-ftT, 0) is very close to the result obtained for the triangular lattice with 
n.n. hopping only.Eil) As for the charge susceptibility, a divergence can be caused only by an attractive Fourier 
component of the potential Vq. Vq has a minimum at the BZ corners ±K, with V±k = — 3V for the n.n. model 
(A) (V±k ~ —1.5422V if a Coulomb tail is added, model B). This minimum leads to an instability towards a 3 x 3 
CDW as {U + 2Vk)Xo(K) = -1, i.e., given our value of % D (K) « 0.2/ti, (U + 2Vk) ~ -5ti. For model A we get a 
transition, when U = 0, at V C MF /t 1 w 0.83. 

In general, the small coupling paramagnetic metal is surrounded by an intermediate coupling region, where com- 
plicated incommensurate - generally metallic - solutions occur. For stronger U and V, commensurate solutions are 
privileged In view of the fact that a 3 x 3 CDW is experimentally relevant, we concentrate our analysis on the 
simplest commensurate phases. These are easy to study with a standard Hartree-Fock (HF) mean-field theory. In 
particular, we restrict ourselves to order parameters associated with non-vanishing momentum space averages of the 
type (cjj ^Ck.o-') and (c£ Ck±K,o-')- Possible non-vanishing order parameters are the uniform magnetization density 
m, 

BZ 



m = jr H( c L(o%/3 c k,/3> = Tf (Stot) , (8) 

S k a,/3 S 

the K-component of the charge density, 

j BZ 

PK = W H( C k,<x C k-K, ff ) , (9) 
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and the K-component of the spin density 



BZ 



'K 



k a,fi 



(10) 



Note that only pk and Sk are 3x3 periodic. Moreover, X-components of bond order parameters of the type 
(cJ a c T ^ a >) are automatically included in the calculation, pk and Sk have phase freedom, and are generally complex: 
Pk = |pK|e l< ^ p , etc. The role of the phase is clarified by looking at the real-space distribution within the 3x3 unit 
cell. For the charge, for instance, (n rj ) = 1 + 2|/5k| cos (2npj/3 + P ), where pj = 0, 1, 2, respectively, on sublattice A, 
B, and C. The e-ph coupling is included but, after linearization, the displacement order parameter is not independent, 
and is given by (zk) = {g/M^vdP'K- Only the phonon modes at ±K couple directly to the CDW. The phonon part 
of the Hamiltonian can be diagonalized by displacing the oscillators at ±K. This gives just an extra term in the 
electronic HF Hamiltonian of the form A?7(/9 k/ ok + H.c), with an energy AU = — g 2 /Mw K which is the relevant 
coupling parameter. This term acts, effectively, as a negative-C/ contribution acting only on the charge part of the 
electronic Hamiltonian. 

With the previous choice of non-vanishing momentum space averages, the Hartree-Fock Hamiltonian reads: 

BZ 



H H -f = ^^ek«-k, CT - Um-S 

k cr 

BZ 

+EE 
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Ck, CT C k+K,<T + H.C. 



U J2{ S K C l,l Ck + K ^ + 5 K C k,T C k+KU + H.C.} 
k 
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EEK CT) 4, CT c k , CT + [4^ c t iCTCk+K) .+H.c.]} 



k cr 
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EE{v T 'wv. 
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(11) 



The last two terms originate exchange contributions due to the V-term; ^ and B^* are shorthands for the 
following convolutions: 



.(cr'cr) _ 



- BZ 



B 



k N. 

1 



k' 

BZ 



s k' 



Vk-k' (Ck'+K,<r C k', CT ' 



(12) 



The BZ is divided into three regions: a reduced BZ (RBZ), and the two zones obtained by k ± K with k 6 RBZ. 
The HF problem in Eq. 11 reduces to the self-consistent diagonalization of a 6 x 6 (including the spin) matrix for 
each k e RBZ. 



A. Landau theory 



The mean-field solutions must be compatible with the symmetry of the problem. Before discussing the HF phase 
diagram we obtain, it is useful to present a few general phenomenological considerations based on a symmetry analysis 
of the Landau theory built from the CDW order parameter /5k (a complex scalar), the SDW order parameter Sk 
(a complex vector), and the uniform magnetization m (a real vector) .E3 In the absence of spin-orbit coupling, the 
possible contributions to the Laundau free energy F allowed by symmetry, up to fourth order, have the form 

F = la p \p K \ 2 + \a m \xa\ 2 + ^a s |S K | 2 + F 3 + F 4 
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^3 = (B p p 3 K + c.c.) + [B ps p K {S K -S K ) + c.c] 

^4 = fe P |p K | 4 + 6 m |m| 4 + 6«|S K | 4 + ^ 2) (S K x S K ) 2 

+ VsIpk| 2 |S k | 2 + & pm |PK| 2 |m| 2 + 6«|m| 2 |S K | 2 + ^(m • S K )(m • S K ) 

+ [b£l(m ■ S K )(S K • S K ) + cc] + [b pmsP 2 K (m • S K ) + c.c] , (13) 

with |Sk| 2 = (Sk • S K ). Notice that third order invariants are present xliie to commensurability, 3K = G (reciprocal 
lattice vector). Therefore, first order transitions are generally possible£3 

This expansion suggests a number of additional comments: i) A CDW can occur without accompanying magnetism, 
i.e., pk 0, while m = and Sk = 0. This is the case, as we shall see later, for the small U region of the HE 
phase diagram, ii) The possible SDW phases are either collinear (1-SDW) (for which (Sk x S k ) = 0) or coplanar.Eil 
The latter have, with a suitable choice of the phases, S K = |Sk| cosa and S K = — i|Sk| sina, and can be generally 
described as a spiral SDW (s-SDW) 

(S r ) = 2|S K | [xcosq cos (K • r) - y sina sin (K • r)] , (14) 

with an eccentricity parameter a ^ 0, 7r/2. (a = or ir/2 are actually 1-SDW along the x or y directions.) a = ir/4 
describes a circular spiral SDW. Now, the only possibility of having a SDW without CDW is via a circular spiral 
SDW (a = 7r/4). Indeed, the third order invariant [-B ps Pk(Sk ■ Sk) + c.c] vanishes by symmetry only for a circular 
spiral SDW, for which (Sk ■ Sk) = 0; in all other cases, a SDW implies - if B ps ^ - a CDW as well, Hi) 
The simultaneous presence of a SDW and a CDW implies, generally, a finite magnetization m, via the fourth order 
invariant [6 pms p K (m • Sk) + c.c], unless the phases of p and S are such that 2cf) p + tt/2 + nn. This happens 

in phase E of our phase diagram, which has therefore no uniform magnetization, iv) The presence of a SDW leads, 

generally, to a finite uniform magnetization as well, via the fourth order invariant [bin S (m. - Sk)(Sk • Sk) + c.c], unless 
the phase 4> a is such that i<p a = n/2 + mir. 



B. Phase diagram in the Hartree-Fock approximation 



We present a brief summary of the mean-field HF calculations for arbitrary U, V, and g, obtained by solving 
numerically the self-consistent problem in Eqs. p|-|l2]. The main phases present in the HF phase diagram are shown in 
Fig. |J for the case of g = 0. The effect of g ^ will be discussed further below. 

Phase A: Spiral SDW insulating phase. The circular spiral SDW (phase A) dominates th e larg e U, small V 
part of the phase diagram, as expected from the Heisenberg model mapping at U — > oo (see sect. Ill A). This is the 
Mott insulator phase, probably relevant for SiC. Its HF bands are shown in Fig. ||(a). 

Phase A': Collinear SDW with m z = 1/3 insulating phase. This is another solution of the HF equations 
in the large U, small V region. It is an insulating state characterized by a linear 1-SDW plus a small CDW with 
4> p = 0, accompanied by a magnetization m z = 1/3 (phase A'). This collinear state lies above the s-SDW by only a 
small energy difference (of order 0.03ii per site), and could be stabilized by other factors (e.g., spin-orbit). A recent 
LSDA calculation for \/3-Si/Si(lll) has indicated this 1-SDW as the ground state, at least if spins are forced to be 
collinear Jl3 The HF bands for this solution are shown in Fig. 0(b), and are very similar to the LSDA surface band for 
Si / Si(l 11). The phase <fi p = of the CDW order parameter corresponds to a real-space charge distribution in which 
one sublattice has a charge 1 + 2|pk|, while the remaining two are equivalent and have charges 1 — \pn\, compatible 
with the experimental findings on Sn/Ge(lll) and Pb/Ge(lll). The amplitude |pk| of the CDW is in general quite 



small i 
density. 



i phase. It should be noted, however, that a STM map is not simply a direct measure of the total charge 
This will be discussed in sect. [v|. 
Phase B': Asymmetric CDW with m z = 1/3 insulating phase. By increasing the n.n. repulsion V, the 
energies of the s-SDW and of the 1-SDW tend to approach, until they cross at a critical value V c of V. At U jt\ = 10 
we find Vc/ti w 3.3 for model A, V c /t\ w 6.6 for model B. As V > V c , however, an insulating asym metric CDW 
(a-CDW) prevails. This is simply the spin collinear version of the non-collinear phase described in Sect. [II B . Fig. || 
shows the energy per site of the most relevant HF solutions at U jt\ = 10 as a function of V for model B (Coulomb 
tail case). The s-SDW and the 1-SDW cross at V c « 6.6£i where, however, the a-CDW insulating solution starts to be 
the favored one. This large- V solution has a large CDW orde r parameter with <f> p ^ (mod. 27r/3), a concomitant 
1-SDW, and m z = 1/3. By recalling the discussion in sect. HI B , we notice that a state with a magnetization m z = 1/3 
and a 1-SDW is the best HF solution once a 3 x 3 restriction has been applied, since a spiral SDW on the singly 
occupied sublattice would involve a larger periodicity (phase B). 

Phase D: Symmetric non-magnetic CDW metallic phase. For small values of U and V, or for large enough 
e-ph coupling g, a metallic CDW with (f) p = (m-CDW) is found. (See Fig. ^(c) for the HF bands.) This phase 
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constitutes a candidate, alternative to the magnetic phase B', and compatible with the main experimental facts, which 
might be relevant for the case of Pb/Ge(lll) and of Sn/Ge(lll). The degree of metallicity of this phase is much 
reduced relative to the undistorted surface (pseudo-gap). 

We stress that the e-ph interaction can stabilize the <j) p = m-CDW also at relatively large U, by countering U with 
a large negative AU = —g 2 /Moj^_. We demonstrate this in Fig. |7|, where we plot the energy per site as a function 
of AU at U/ti = 8 and V/tx = 2, for the three relevant HF solutions, i.e., the spiral SDW (phase A), the collinear 
SDW with m z = 1/3 (phase A'), and the metallic non-magnetic CDW (phase D). The spiral SDW is unaffected by 
the electron-phonon coupling. The energy of the collinear SDW with m z = 1/3 improves a little bit by increasing 
g, due to the small CDW amplitude of this phase. This effect is not large enough as to make this phase stable in 
any range of couplings. At a critical value of g, the metallic non-magnetic CDW (where the CDW order parameter 
is large, |/?k| ~ 0.5) wins over the magnetic phases. The Fourier transform of the lattice distortion at K is given by 
(z K ) = {g/Muj Si )p K =p b f K \AU\/g. 

A rough estimate shows that the order of magnitude of the electron-phonon coupling necessary to stabilize the 
CDW phase is not unreasonable. With g — 1 eV/A, Msi — 28, and o>k ~ 30 meV we get AU w — 3ti, sufficient to 
switch from a s-SDW ground state to a m-CDW for U/t± = 8 and V/t± = 2. With these values of the parameters 
we have |pk| ~ 0.43, and we estimate |(zk)| ~ 0.07A. This corresponds, since (z r ) ~ 2cos(K • r)|(zK)|, to a total 
displacement between the adatom going up and the two going down of Az = 3|(zk)| ~ 0.2A. 

We notice that values of g much larger than those used in Fig. ^ would eventually stabilize a superconducting 
ground state (see Appendix). 



V. CDW ORDER PARAMETER AND STM EXPERIMENTS 

We discuss, in the present section, the relationship between the CDW order parameter, as defined in Eq. ^, and an 
STM map of the surface. As the crudest approximation to the tunneling current for a given bias Vbias we consider the 
integral of tip, charge density for one-electron states within Vbias from the Fermi level, weighted with barrier tunneling 
factor T{V)fyB 

/•V biaB 

J{V hias ,r = x,y;z)& dV V |*„ k (r)| 2 <5(^„ k - E F + V)T(V) . (15) 

Jo nk 

The tunneling factor leads to weighting prominently the states immediately close to the Fermi level. In view of the 
purely qualitative value of Eq. (|l^), we have moreover decided to ignore T(V) altogether and to account for its effect 
by reducing the bias voltage V e s in Eq. (|l5|), to an effective value Iq^as- By doing this, we have extracted an "STM 
map" for a point in phase A' (U/ti — 9 and V = 2, model A) - a spin-density waves where the amplitude of the CDW 
order parameter is rather small, |pk| = 0.039 - and a point in phase D (U/ti = 4 and V = 2, model A) - a pure CDW 
where the order parameter is quite large, \pn\ — 0.4. The results for constant z, and x, y moving from adatom A to 
B to C, are shown in Fig. ||(a) and (b), for the two cases. The solid curves refer to positive bias (current flowing from 
the sample to the tip), probing occupied states close to the Fermi level. The dashed curve refers to negative bias, 
probing unoccupied states. In both cases a) and b), one of the three atoms yields a larger current at positive bias, 
while the other two atoms have larger currents at negative bias. The insets show the predicted "contrast" between 
the two peak values, (Ji — J2)/(Ji + J%), J\ and J 2 being in each case, respectively, the largest and the smallest of 
the STM peak currents at the positions of the adatoms. We notice the following points: i) for the occupied states 
(positive bias) the pure CDW phase has, as expected, a larger contrast than the magnetic phase. As we neglect the 
tunneling factors T(V), in the limit of large positive effective bias we recover the total asymmetry in the charge of 
the two inequivalent atoms, (ni — n-i)j(n\ + n^), indicated by a dashed horizontal line in the insets. Observe that the 
way this large bias limit is reached is completely different for the two cases a) and b): in the magnetic case a) the 
contrast overshoots at small biases attaining values substantially larger than the nominal CDW order parameter, and 
then goes to the limit (n\ — n^)/ (rii + n^j from above; in the pure CDW case b), on the contrary, the limit is reached 
monotonically from below, ii) for empty states (negative bias) the contrast is even more surprising: at small bias it is 
very large in both cases a) and b). By increasing the bias, the contrast for the pure CDW case tends monotonically 
to a large value, whereas the magnetic case shows a strong non monotonicity. 

These results suggest that one should look more carefully, and quantitatively, at the behavior of the asymmetry 
between STM peak currents as a function of the bias, including the region of relatively small biases: the different 
behavior of the asymmetry of the magnetic case versus the pure CDW case should be marked enough - and survive 
in a more refined analysis including T(V) - as to make the STM map a good way of discriminating between the two 
scenarios. 
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VI. DISCUSSION AND CONCLUSIONS 



Within our model study we have learned that on the surfaces considered: 

(i) If U and V are ignored, there is no straight electron-phonon driven 3x3 CDW surface instability. However, any 
phase involving a CDW, for example as a secondary order parameter attached to a primary SDW, can take advantage 
and gain some extra stabilization energy from a small surface lattice distortion, via electron-phonon coupling. 

(ii) Electron-electron repulsion and the two-dimensional Fermi Surface are capable of driving transitions of the 
undistorted metallic surface to a variety of states, that are either insulating or in any case less metallic, some possessing 
the 3x3 periodicity. 

(iii) This can occur via two different mechanisms: a) the inter-site repulsion V can stabilize insulating or semi- 
metallic CDWs, without a crucial involvement of spin degrees of freedom; b) the on-site repulsion U can produce 
essentially magnetic insulators with or without a weak accompanying 3x3 CDW, as required by symmetry. 

(iv) For U moderate of order W and for smaller V, an interesting state is realized, with a large SDW and a small 
accompanying CDW. The state is either a small-gap insulator, or a semi-metal, and may or may not be associated 
with a net overall magnetization, depending on the nature (linear or spiral, respectively) of the leading SDW. 

(v) For U and V both small but finite, a metallic CDW without any magnetism is obtained. The same phase can 
also be stabilized for larger values of U by the presence of a substantial electron-phonon coupling. We stress that, in 
this case, V is the coupling responsible for the 3x3 symmetry of the unit cell, whereas the role of the electron-phonon 
coupling is that of destroying magnetism by effectively decreasing U. Electron-phonon coupling alone is not sufficient 
to justify a commensurate 3x3 CDW. 

(vi) Either of the phases in (iv) or (v) could be natural candidates for explaining the weak 3x3 CDW seen 
experimentally on Sn-Pb/Ge(lll). 

(vii) Finally, for large U, small V (in comparison with the bandwidth W) the Mott-Hubbard state prevails. It is a 
wide-gap insulator, with a pure spiral SDW, with 3x3 overall periodicity, and coplanar 120° long-range spin ordering 
at zero temperature. It possesses no net magnetization, and no accompanying CDW. 

(viii) The above is the kind of state which we expect to be realized on SiC(OOOl), and also possibly on K/Si(lll):B. 
Among existing experiments, we have addressed particularly photoemission and STM. Our calculated band structure 

for both the SDW/CDW state A' (iv) and th&jpure CDW state D (v) exhibit features which are similar to those 
found in photoemission from Sn-Pb/Ge(lll).B~En3 The simulated STM images for the two kind of states are predicted 
to differ in their voltage dependence. 

Future experiments are strongly called for, aimed at detecting whether magnetic correlations are actually dominant, 
as we think is very likely, on all these surfaces, or whether Sn-Pb/Gc(lll) are instead non-magnetic and electron- 
phonon driven. The issue of whether magnetic long-range order - which we definitely propose for SiC(OOOI) and 
K/Si(lll):B at T = 0, and also hypothesize for Sn-Pb/Ge(lll) - survives up to finite temperatures is one which we 
cannot settle at this moment. This due to the difficulty in estimating the surface magnetic anisotropy, without which 
order would of course be washed out by temperature. In any case, it should be possible to pursue the possibility of 
either magnetism or incipient magnetism using the appropriate spectroscopic tools. 

This line of experimental research, although undoubtedly difficult, should be very exciting since it might lead to 
the unprecedented discovery of magnetic states at surfaces possessing no transition metal ions of any kind, such as 
these seemingly innocent semiconductor surfaces. 

We acknowledge financial support from INFM, through projects LOTUS and HTSC, and from EU, through ER- 
BCHRXCT940438. We thank S. Modesti, J. Lorenzana, M.C. Asensio, J. Avila, G. Le Lay, EW. Plummer and his 
collaborators, for discussions. 



VII. APPENDIX. LARGE NEGATIVE U: A SUPERCONDUCTING GROUND STATE. 

The limit of large negative U, U — > — oo, is considered hepe to show that CDWs are not favored by on-site attraction 
alone. Instead, a superconducting ground state is favored.Ej To see this, consider the real-space states which are the 
low energy configurations for U — > — oo: they consist of N e /2 sites (if N e is the number of electrons) each of which is 
occupied by a pair of electrons with opposite spins. The large degeneracy in this manifold of states is - once again, like 
in the U — > oo case - removed by kinetic energy in second order perturbation theory. By assigning a pseudo=spin-l/2 
state to each site (up, if occupied by a pair, down if empty) one can show that the effective Hamiltonian iai3 

^-Ey ( s rA~ + H - c -) + E J ij - ( 16 ) 
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with Jjj = 4\tij\ 2 /\U\ and Jfj = J^. If T^-terms are added, J z is modified to Jfj = + 4Vij. Restricting our 
consideration to the n.n. case, we are left with a n.n. Heisenberg Hamiltonian with ferromagnetic xy-part and an 
antiferromagnetic z-part. The sign of the xy-part cannot be changed at will by a canonical transformation because 
the lattice is non-bipartite. The result is that the order is in the plane (i.e., superconductivity wins) for small V. 
Only if V is large enough the CDW (i.e., order in the z-direction) will be favored. 

Entirely similar considerations apply to the case of strong electron-phonon coupling, g — > oo. 
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Figure Captions 



FIG. 1. Surface state dispersion for hypothetical Si/Si(lll), as obtained from gradient corrected LDA (solid squares). A 
very similar band is obtained for Pb/Ge(lll) and Sn/Ge(lll). The solid line is a tight-binding fit obtained by including up 
to the sixth shell of neighbors, ti,--- ,t 6 . The fit gives ty = 0.0551 eV, and t 2 /ti = -0.3494, t 3 /ti = 0.1335, U/ti = -0.0615, 
ts/ti — 0.0042, ts/ti — —0.0215. The dashed line is the best fit using ti and ti only. Upper inset: The Fermi surface of the 
half-filled surface band. The outer hexagon is the BZ of the y/3 x y/3 phase, and the inner hexagon is the BZ of the 3x3 phase. 
Notice the poor nesting at the BZ corner wavevector K = (47r/3a, 0), joining two opposite Ms X 3 points. Lower inset: The zero 
temperature Lindhard response function x°(q) f° r the half-filled surface band. Notice the two peaks located at qi « 0.525K 
and q2 fa 1.32K, and no feature whatsoever at K, indicating poor nesting. 

FIG. 2. Density contours of the Wannier function associated with the Si/Si(lll) surface band, calculated with gradi- 
ent-corrected LDA: dots correspond to atomic positions. 

FIG. 3. The CDW on the adatom triangular lattice in the limit fey = 0, U <C V . Sublattice A is doubly occupied, sublattice 
B is singly occupied, and C is empty. The large spin degeneracy associated to the unpaired singly occupied sites on sublattice 
B is removed by the next-nearest neighbor hopping ti in second order perturbation theory. The arrow indicates a possible 
virtual process, leading to an extra doubly occupied site (with an associated energy gap of U), which generates the standard 
antiferromagnetic exchange coupling between the spins on sublattice B. 

FIG. 4. Schematic Hartree-Fock phase diagram of model A (nearest neighbor V only, g = 0) for the band structure shown 
in Fig. [lj Only the most important commensurate 3x3 phases have been studied. (Details of the merging of the various lines 
are not accurate.) Left figure: phase obtained allowing spin non-collinearity. Right figure: strictly collinear phase diagram. 
The non-vanishing order parameters of the different phases are as follows: A: Sk = |Sk|/v2(1, —i, 0); A': Sk = |Sk|(0, 0, 1), 
m z = 1/3, p K = |pk|; B ! : S k = |S K |(0,0, 1), m z = 1/3, p K = IpK^" (with < <j> p < tt/6); B: 3^ extension of phase 
B' when allowing for non-collinearity of the unpaired spins. C and C are semi-metallic versions of A and A', respectively. 
D: pK = |px|; E: Sk = i|Sk|(0, 0, 1), pn = — |pk|; F: pk = |pK|e l< ^ p (with < 4> p < n/6); IM: Incommensurate metallic 
SDW/CDW; PM: Paramagnetic metal. Phases A(A'), and B(B') are insulating and magnetic. Phases C(C'), E, and IM are 
metallic and magnetic. Phase D is a pure CDW and is metallic. Phases A', C', and D have CDW order parameter with 
the same symmetry as observed on Pb/Ge(lll) and Sn/Ge(lll). The effect of a finite electron-phonon coupling (g 7^ 0) is 
discussed in the text. 

FIG. 5. Plot of the HF electronic bands along high symmetry directions of the BZ for the s-SDW and two CDW <j> p — 
solutions: (a) at U/ti = 9 and V/ti = 2, the insulating s-SDW (phase A, ground state); (b) at U/ti = 9 and V/ti = 2, the 
insulating solution with a small CDW and m z = 1/3 (phase A', meta-stable, the actual ground state being the s-SDW); Solid 
and dashed lines denote up and down bands, respectively, (c) at U/ti = 4 and V/ti = 2, the metallic solution with a large 
CDW and no magnetism (phase D). The band structure for phases C and C are similar to A and A' except for band overlap, 
making them semi-metallic. Insets indicate the charge and spin imbalance (when present) between the three adatoms in the 
3x3 unit cell. 



FIG. 6. Energy per particle, as a function of V at U/tj = 10, for the commensurate HF solutions of phase A, A', and B', 
obtained for model B (Coulomb tail case). The results for model A are entirely similar, with V c w 3.3ii. 

FIG. 7. The energy per site, as a function of the electron-phonon coupling — At/ = g 2 /Mu|, at U = 8ti and V — 2t\ for the 
model with n.n. interactions only, for the three relevant HF solutions, i.e., the spiral SDW, the collinear SDW with m z = 1/3, 
and the metallic non-magnetic CDW. 

FIG. 8. (a) STM map for tip motion (at constant height) along the triangle indicated, for the 1-SDW state at U/ti = 9, 
V = 2. The result is obtained from Eq. [IH] without explicitly including T(V), but reducing the bias to an effective one, 
^bSs — 6ti, thus imposing that only states sufficiently close to the Fermi level contribute significantly to the STM current. The 
solid and dashed curves refer, respectively, to positive and negative STM bias. The inset shows the "contrast" between the 
peaks (see text) as a function of the effective bias V^ B . (b) Same as (a) for a pure CDW state in phase D. 
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